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We study the prospect of dark matter (DM) searches in the monojet channel at 
future pp colliders with center-of-mass energies of 33, 50, and 100 TeV. We con¬ 
sider a class of simplified models in which a vector boson connecting DM particles 
to quarks is introduced. Comparing with studies in the effective held theory, the 
present framework gives more reasonable production rates and kinematics of the 
DM signatures. We estimate the sensitivities of future colliders with an integrated 
luminosity of 3 ab^ 1 to the DM-induced monojet signature and show the parameter 
space that can be explored. The constraints from direct and indirect DM detection 
experiments are compared with the future collider sensitivities. We find that the 
future collider detection will be much more sensitive than the indirect detection for 
the vector interaction, and have better sensitivities than those of the direct detection 
by several orders of magnitude for the axial vector interaction. 

PACS numbers: 95.35.+d, 12.60.-i 

I. INTRODUCTION 

Observations of astrophysics and cosmology reveal that the main component of matter 
in the Universe is dark matter. However, the nature of DM particles is still unclear. Al¬ 
though the standard model (SM) has achieved great success after the LHC discovery of a 
~ 125 GeV Higss boson, it cannot provide any suitable candidate for the cold DM. There¬ 
fore, the existence of DM gives a hint of new physics beyond the SM (BSM). The most 
popular and attractive candidates for DM are the so-called weakly interacting massive par¬ 
ticles (WIMPs), which can be thermally produced in the early Universe and naturally give 
the correct observed DM relic density. Because of their weak interactions with SM particles, 
WIMPs with masses of ~ (9(10 2 ) GeV are expected to be produced at high energy collid- 
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ers. Therefore, searching for such DM particles is a very important task for future collider 
experiments. 

DM particles can be produced by cascade decays of heavy new particles in new physics 
models at colliders. A typical example is the lightest neutralino in supersymmetric models, 
which predict a bunch of new particles that have not been discovered by LffC searches up to 
now. Another possibility is that DM particles are directly pair-produced in collisions of SM 
particles. In this case, DM signatures have been widely studied in the context of effective held 
theory (EFT) 1-10], which provides a “model-independent” way to quantitatively compare 
the sensitivity among collider, direct, and indirect detection experiments. Since DM particles 
escape from the detector without energy deposit, an additional energetic jet or photon 
is required to reconstruct the signature with missing transverse energy (^ T ). The latest 
ATLAS and CMS searches have not found such so-called “monojet” and “monophoton” 
signatures induced by DM particles, and hence set constraints on the energy scales of effective 
contact operators describing DM interactions with SM particles in the EFT llNl4|. 

As long as the mediators interacting with both DM and SM particles are so heavy that 
they can be integrated out, the EFT approach is valid. However, it has many shortages 


and would break down in many cases, as pointed out in several recent works 
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23]. It is 


well known that the EFT fails when the typical momentum transfer involved in the reaction 
is comparable to the mediator mass. This means that in order to safely use the EFT, the 
scale of BSM physics should be much higher than the collision energy, otherwise mediators 
may be directly produced in collisions. Besides, the EFT approach is invalid for a mediator 
with a width comparable to or larger than its mass la, [23]. Furthermore, unitary and 
perturbativity conditions could also set constraints on the validity of the EFT fril |. 

In a more appropriate approach, so-called “simplified models”, mediators with moderate 
masses are introduced to connect DM particles to SM particles. In principle, the simplified 
model approach can be mapped into the EFT approach, as studied in the limit of heavy 
mediators 5, 7, 9]. On the other hand, it can be realized as particular cases of UV 
complete models, and has been widely used in supersymmetry studies (see e.g. Refs. Q- 


26]). For a mediator with a moderate mass, full kinematics and topologies of DM signatures 


at colliders can be studied in details. Furthermore, collider constraints and reaches could 


also be easily compared with 


simplified models 18l. 1231,127-134] 


hose from direct and indirect detection experiments in specific 
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In this work, we study a class of minimal simplified models involving a DM particle and 
a neutral spin-1 mediator Z'. Each model is characterized by the following parameters: the 
mediator mass my, the DM particle mass m x , the Z' coupling to the DM particle g x , and 
the Z' couplings to quarks g q . Recent 7 and 8 TeV LHC results can be used to set limits on 
the parameter space of these simplified models 31]. It is expected that LHC searches 
will set more stringent constraints with larger integrated luminosity at its designed energy 
13 — 14 TeV in the next few years. In this work, we investigate the prospects of exploring 
DM signals in the framework of simplified models with Z' at future pp colliders. Three 
collision energies of 33, 50, and 100 TeV for the very large hadron collider and the super 
proton-proton collider are considered (for recent phenomenology studies, see Refs. 35-139]). 


We study the sensitivity of future colliders to the monojet signatures and set constraints on 
the parameter space. It is expected that these future colliders will be powerful to detect DM 
particles and mediators with masses of (9(TeV). Finally, we compare the reaches of future 
collider searches with results from direct and indirect detection experiments. 

This paper is organized as follows. In Sec. II, we describe the minimal simplified DM 
models containing a DM particle and a vector mediator. In Sec. Ill, we discuss the kinematics 
of the monojet signature and estimate the sensitivity to these models at future colliders. In 
Sec. IV, we present constraints and future reaches of direct and indirect searches, and then 
compare them with the capability of future collider. Finally, we conclude in Sec. V. 


II. MINIMAL SIMPLIFIED DARK MATTER MODELS 

In this section, we describe a class of minimal simplified models, in which a spin-1 media¬ 
tor Z' connecting the DM particle (y) to quarks is introduced. The DM particle is assumed 
to be a singlet under SM gauge symmetries and only couples to Z'. Here we are not going 
to enumerate all the possible sets of particle spins and Lorentz structures; we only focus on 
two typical cases where the DM particle is either a Dirac fermion or a complex scalar. For 
the fermionic DM case, we consider DM particles and quarks have a vector or axial vector 
interaction with Z\ which reads 

£ fv = ggZ'^q + g x zlxfx, ( 1 ) 

£ F a = 9qZ'^qYi b q + g x z '^xYi b Xi 

q 


( 2 ) 
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where the sums run over all quark flavors. For the scalar DM case, axial vector current cannot 
be constructed for the DM particle, thus we only consider the vector current interaction given 
by 

£sv = dqZ'^qYq + igxZ'tfoC^X ~ ^X*)x\- (3) 

q 

The subscripts FV and FA represent the fermionic DM case with vector and axial vector 
interactions, respectively, while SV represents the scalar DM case with vector interaction. 
Here we have assumed that all the interactions are renormalizable and CP invariant. These 
simplified models are typical, and their discussions can be easily extended to other models 
with other interactions in principle. 

Free parameters in these models are m x , m Z r , g X i and g q ■ In general, the Z' couplings to 
quarks g q are flavor dependent, and should be determined by the details of underlying UV 
complete models. For simplicity, we just assume g q is universal for all the quarks; such an 
assumption can provide an illuminating example of DM searches. 

Comparing with g x , g q would suffer from additional limits from dijet resonance searches 
for qq —>• Z' —)■ qq at hadron colliders 40|, l4l|. Especially, it is expected that future colliders 
will have strong capability to search for dijet resonances due to their high energies and 
luminosities. Combining the results from dijet and monojet searches, constraints on g q and 
g x can be simultaneously set. In this work, we only focus on monojet searches and simply set 
g q = g x in order to further reduce the number of free parameters. Since the DM production 
rate is almost independent from the Z' width Y z > and just sensitive to the product of g q 
and g x in most of the parameter space except for the resonance region, the discussions on 
monojet searches under the assumption g q = g x can be easily extended to the g q ^ g x case 
in these regions. 

r z , would play an important role in the resonance region. Since Y z > depends on the 
sum of the terms proportional to g q and g x , Y z > and g q g x can be simultaneously treated as 
free parameters instead of g q and g x . The detailed discussions on the effect of Y z > and the 
g q 7^ g x case in DM searches can be found in several works [16, Q, 
can be expressed as 


M 


42]. The Z' width 


F z> = r {Z' ->• XX/xX*)®(rn Z ' - 2 m x ) + ^ c q Y(Z' ->• qq)Q(m z > - 2m,), (4) 

where the step function means that the particular Z' decay channel opens when it is allowed 
by kinematics. The color factor c q = 3. Partial widths contributed by individual decay 
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channels are 


r F v(Z' —* qq) 
r F A (Z' —> qq) 
rs v{Z' —> qq) 


™z’ 2 C (i . 
12n 9q ^ q \ + 

m Z' 2 r 3 

12n 9 ™’ 
r F y(Z' —> qq), 



r FV (z' —* xx) 
r FA (z' ->• xx) 
r sv(z’ xx*) 


m z > 2 
127 t 9x 


9\^x 


m z > 

12 vr 


q 2 C- 


2 ITT 

1 H-— 

mi, 


7 J}tz_9xc 3 

12vr 4 


(5) 

( 6 ) 
(7) 


Here £/ = y 1 — 4mj/m 2 z , and can be regarded as the particle velocity in the Z' rest frame. 

If the Z' width is larger than the Z' mass, it would be questionable to treat Z' as a 
particle. Therefore, the couplings g x and g q should not be too large. In view of the fact that 
< 1 and £ g (l + 2m 2 q /m 2 z/ ) < 1, we have Yz> < {g x + J2 q c q gq) mz'/{Y2n) for the fermionic 
DM case. Consequently, the requirement Yz> < gives 


9x + c i g g < 127T ■ ^ 8 ) 

q 

Under the assumption of g x = g q , we can get a rough limit of g q < 1.4. The limit on g q for 
the scalar DM case is similar. Moreover, the perturbativity condition also requires that the 
couplings should not be too large. For these reasons, below we only consider the parameter 
space with g q (g x ) < 1. 


III. MONOJET SIGNATURE AT FUTURE COLLIDERS 


If DM particles are produced in pair via high energy pp collisions, an additional energetic 
object, such as a jet, photon, charged lepton, or Z boson, is required to reconstruct |5 T 
and trigger the event. In the following, we investigate the monojet signature pp —» > 

XX/XX*) + jets at future pp colliders with center-of-mass energies of 33, 50, and 100 TeV. 
Background and signal samples at the parton level are generated by MadGraph 5 43], to 
which the simplified models are added through FeynRules 2 44], We use PYTHIA 6 45] 
to deal with particle decay, parton shower, and hadronization processes, and adopt the 
MLM matching scheme. Delphes 3 46( is utilized to carry out a fast detector simulation. 
Jets are reconstructed using the antialgorithm 47] with a distance parameter R = 0.4. 
For realistic search strategies using a future detector designed with higher efficiency and 
resolution, results in this work are expected to be improved. 

The dominant SM backgrounds are Z(—» uu) + jets and W{—> Iv) + jets. For the W{—> 
Zz/)+jets process, charged leptons may be clustered into a nearby jet or undetected when they 
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close to the beam pipe. On the other hand, the Z(—> m/)+jets background is irreducible. We 
optimize selection criteria to efficiently suppress backgrounds and maximize the statistical 
significance. 

There should be at least an energetic jet in the final states. The leading jet ji is required 
to have 1 77 (jT)| < 2.4 and Pt(ji) > 1.6/1.8/2. 6 TeV for y/s = 33/50/100 TeV. Events with 
more than two jets with p t > lOOGeV and \rj\ < 4 are rejected. The DM production process 
may involve more than one jet from initial state radiation. In order to select more signal 
events, a second jet ( 72 ) is allowed if it satisfies the condition A<f>(ji, j 2 ) < 2.5. The cut on 
A0(ji, j r i) is necessary to suppress the QCD multijet background, where large fake may 
come from inefficient measurement of one of the jets. Furthermore, in order to reduce other 
backgrounds, such as W{—> Iv ) + jets, Z(—> W _ ) + jets, and tt + jets with leptonic top 
decays, the events containing isolated electrons, muons, taus, or photons with p T > 20 GeV 
and |?/| < 2.5 are discarded. 

The cut condition on is crucial for DM searches. I 11 Fig. [lj we show the normalized 
distributions of signals and backgrounds, assuming g q = g x = 0.1, m x = 1 TeV and 
mz' = 5 TeV. The distributions become harder as the collision energy increases. The slopes 
of the distributions are mainly determined by the energy scales involved in collisions. Since 
here mz' and m x are chosen to be much larger than mz and rriw, signal distributions are 
always harder than background distributions. The distributions for three simplified models 
are similar, implying that spins and Lorentz structures play minor roles in high energy 
collisions. I 11 order to optimize the statistical significance, we choose the cut condition as 
> 1.6/1.8/2.6 TeV for = 33/50/100 TeV. 

TABLE I: Cross sections (in fb) after cuts for backgrounds and signals. For DM production signals 
pp —» XX./XX*) + jets, g q = g x = 0.1, m x = 1 TeV, and mz' = 5 TeV are assumed for all 

the three simplified models. 



W(—>• Iv) + jets 

Z(—t vv) + jets 

FV 

FA 

SV 

33 TeV 

8.179 x 10 1 

1.948 x 10 2 

3.043 x 10" 2 

2.399 x 10" 2 

6.133 x 10~ 3 

50 TeV 

6.991 x 10 1 

1.816 x 10 2 

9.037 x 10" 2 

7.054 x 10" 2 

1.824 x 10~ 2 

100 TeV 

3.475 x 10 1 

1.062 x 10 2 

2.340 x 10 -1 

1.851 x 10 -1 

4.735 x 10~ 2 
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(a) y/s = 33 TeV 


(b) v/i = 50 TeV 



(c) ^ = 100 TeV 

FIG. 1: Normalized distributions for signal and background samples at y/s = 33 TeV (a), 
50 TeV (b), and 100 TeV (c) with thresholds of 800, 1000, and 1600 GeV, respectively. For DM 
production signals pp —>• XX/XX*) +jets, g q = g x = 0.1, m x = 1 TeV, and mz> = 5 TeV are 

assumed. The dot-dashed vertical lines denote the locations of the cuts used in the analysis. 

In Table [fl we list the cross sections of backgrounds and signals after cuts at y/s = 33, 
50, and 100 TeV, assuming g q = g x = 0.1, m x = 1 TeV, and mz> = 5 TeV. Since we 
choose stricter cuts at higher collision energy, the cross sections of backgrounds decrease as 
y/s increases. DM production rates in the FV and FA models are similar, implying that 
the influence of Lorentz structures is slight at high energy. Cross sections of scalar DM 
particles are smaller than those of fermionic DM particles by a factor smaller than 10. This 
is because scalar particles have less hclicity states and their production cross section suffers 

























a kinematic suppression in the angular distribution. 




(a) Fermionic DM with vector Z' (b) Fermionic DM with axial vector Z' 



(c) Scalar DM with vector Z' 

FIG. 2: Estimated 90% C.L. limits on the FV (a), FA (b), and SV (c) models in the m x -g q plane 
for the monojet + channel at pp collides with y/s = 33 TeV (blue lines), 50 TeV (purple lines), 
and 100 TeV (red lines), assuming an integrated luminosity of 3 ab^ 1 . Solid and dashed lines 
correspond to mz 1 = 1 TeV and 5 TeV, respectively. The regions above the curves are expected 
to be excluded. We have assumed g x = g q . 

Now we investigate the sensitivity to the simplified DM models at future pp colliders. 
The estimated 90% C.L. exclusion limits in the rn x -g q plane for mz> = 1 and 5 TeV are 
shown in Fig. [21 The DM production cross section dramatically depends on whether the 
mediator Z' is on-shell or not. When m x < mz>/ 2, Z' can be on-shell produced and then 
decays into a pair of DM particles. In this case the DM production cross section is resonantly 
enhanced and proportional to g^Br(Z f —> xx) under the narrow width approximation. For 
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mz> 2m x , the branching ratio to DM particles Bi :{Z' —> xx) is almost a constant, and 
can be estimated as 1/19, 1/19, and 1/73 for the FV, FA, and SV models, respectively. In 
this limit, the DM production cross section is irrelevant to rn x . This explains the behaviors 
of the limits for m x < rriz '/2 in Fig. [2j Note that Br [Z' —>■ xx) m the FV model almost 
equals to that in the FA model. This is why the limits for these two models are nearly the 
same. For the SV model, Br (Z 1 —> xx) is smaller than those of fermionic DM models by a 
factor of ~ 4, hence the limits on g q is weaker by a factor of ~ 2. 

The sensitivity drops quickly as m x becomes larger than mz'/ 2 and Z' is off-shell pro¬ 
duced. In this case, the DM production cross section is proportional to [g q g x /(Q 2 — m|,)] 2 , 
where Q 2 is the typical momentum transfer to the DM particle pair. As a result, for a large 
enough m x , the collider search may be more sensitive to a heavier on-shell Z' than a lighter 
off-shell Z', as shown in Fig. [2] 

When Z' is off-shell produced and m|, <C Q 2 , the DM production cross section scales as 
(dqQx/Q 2 ) 2 1 which is irrelevant to rriz'- It can be found in Fig. [2] that when m x increases, 
the limits for different rnz< become close to each other. On the other hand, if m 2 z , Q 2 , 
the cross section is proportional to (g q g x /m 2 z ,) 2 and can be matched to that in the EFT 
approach with an effective energy scale of A eff = m Z ' / ^g q g x - 

IV. COMPARISON AMONG DIFFERENT DM DETECTION EXPERIMENTS 

DM signal features in different simplified models are quite similar at colliders, but may 
be very different in other DM detection experiments. In this section, we will discuss the 
prospects of different kinds of DM searches and give a comparison among them. Firstly, 
we consider DM direct detection experiments, which search for nuclear recoil signatures 
induced by DM-nucleus scatterings. Since the typical DM velocity near the Earth is ~ 10~ 3 , 
the momentum transfer in the DM-nucleus scattering is ~ 0( KeV), which is much lower 
than that at colliders. Thus the nonrelativistic limit is valid and different DM simplified 
models could induce distinct phenomenologies in the direct detection. For instance, vector 
and axial vector DM-quark interactions induce spin-independent (SI) and spin-dependent 
(SD) DM-nucleus scatterings, respectively. It is well known that the SI scattering cross 
section is coherently enhanced by the square of the nucleon number in the nucleus. On the 
other hand, the SD scattering signature has no coherent enhancement and depends on the 
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particular spin property of the target nucleus. 




(a) FV model, spin independent (b) FA model, spin dependent 



(c) SV model, spin independent 

FIG. 3: Estimated 90% C.L. limits on the FV (a), FA (b), and SV (c) models from the monojet 
searches in the m x -a X N plane. Blue/purple/red lines correspond to yfs = 33/50/100 TeV, while 


solid (dashed) lines correspond to mz> = 1 (5) TeV. 
bounds from direct detection experiments XENON100 


constraints from SIMPLE 


or the SI DM-nucleon scatterin g, r ecent 


M, LUX 


49[|, and SuperCDMS 


50], and 


the expected reach of XENON1T [5j| are shown. For the SD DM-nucleon scattering, current 

52], PICASSO 53 J], and COUPP 54] are shown, as well as the limits 




from neutrino detection experiments Super-K [55] and IceCube [5a]. 


In order to compare reaches of the direct detection and the future collider detection, we 
show in Fig. [3] the translated 90% C.L. limits from the monojet searches in the m x -cr X N 
plane, where <j x n is the DM-nucleon scattering cross section. For the vector interaction, 
current direct detection experiments have set stringent constraints on the SI signature. The 
collider detection would have a stronger capability to search for light DM in the region 
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of m x < 10 GeV, where the direct detection dramatically lose the sensitivity due to the 
experimental threshold. For the axial vector interaction, constraints from direct searches 
are very weak. The monojet search will significantly improve current direct detection limits 
by several orders of magnitude. 


pp collider, 3abd , Vector Z’ 



FIG. 4: Estimated 90% C.L. limits on the FV model from the monojet searches in the m x - 
{(J ann f) plane. Blue/purple/red lines correspond to yfs = 33/50/100 TeV, while solid (dashed) 
lines correspond to rriz' = 1 (5) TeV. For a comparison, also shown are the recent limit from the 


Fermi-LAT 4-year gamma-ray observations on the dwarf galaxies 


57] and the expected limits from 


the AMS-02 anti-proton detection for 20 years and the CTA observation on the dwarf galaxy Segue 1 
581 ]. The horizontal line denotes the canonical thermal relic value 3 x 10~ 26 cm 3 /s. 


for 100 hours 


Indirect detection experiments search for high energy cosmic rays, gamma rays, and 
neutrinos induced by DM annihilations in the Galaxy and extragalactic objects. Signa¬ 
tures in these experiments depend on the thermally averaged DM annihilation cross section 
(<7 ann u) and DM density distributions in annihilating regions. In Fig. [I] we demonstrate 
the translated 90% C.L. limits from the monojet searches in the rn x -(cr ann v) plane. The FV 
model leads to s-wave DM annihilations into quarks and can be explored by gamma-ray and 
cosmic-ray observations. For m x < mz'/ 2, future colliders could have a sensitivity better 
than indirect searches by several orders of magnitude. In the FA model, DM annihilations 
into quarks are helicity suppressed in s-wave and hence depend on the DM velocity disper¬ 
sion ( v 2 ), which is typically ~ 0( 10~ 6 ) in the Galaxy. In the SV model, DM annihilations 
into quarks are of p -wave in the leading order and also highly suppressed by (v 2 ). There- 
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fore, these two models cannot be explored via qq channels in indirect detection experiments. 
Nonetheless, when m x > mz 1 -, DM particles can annihilate into Z' pairs, which would decay 
into quarks and give rise to detectable gamma-ray and cosmic-ray signals. 



(a) FV model 


(b) FA model 



(c) SV model 

FIG. 5: Estimated 90% C.L. limits on the FV (a), FA (b), and SV (c) models in the mz'-m x plane 
for the monojet + channel at a pp collider with yfs = 50 TeV, in comparison with limits from 
direct searches. An integrated luminosity of 3 ab^ 1 is assumed. Dashed blue/purple/green/red 
curves correspond to g q = g x = 0.1/0.3/0.5/1.0. The regions below these curves are expected to 
be excluded. Limits from direct searches (solid lines) are converted assuming g q = g x = 0.5. The 
unitarity violation region for g q = g x = 1 in the FA model is indicated by light red color. 

In the above study, we have fixed rriz 1 ■ Now we discuss the influence of the Z’ mass for 
DM searches. Fig. [5] shows the estimated 90% C.L. limits in the mz'-m x plane from the 
monojet search at a future pp collider with fixed values of g q . The collision energy and the 
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integrated luminosity are taken to be 50 TeV and 3 ab -1 , respectively. The parameter space 
can be reasonably divided into two parts, m x < mz'/2 and m x > mz'/ 2. As mentioned 
above, the DM production cross section is controlled by the factor [ g q g x /(Q 2 — m|,)] 2 in the 
region of m x > mz> /2,. For m x mz> and thus Q 2 m|,, the DM production rate would 

be irrelevant to rn z 1 • This explains the nearly horizontal segments of the estimated limits, 
except for the limits in the FA model when mz> < 50 GeV, which will be discussed later. 


For a weak coupling like g q = g x = 0.1, the DM production rate with off-shell Z' may be too 
low to detect, hence the limits are restricted to the region near m x = mz>/ 2. In the region 
of m x < mz'/^i Z' is on-shell produced and the limits are nearly vertical. The maximum 
reach of rn zj only depends on the collision energy and is almost irrelevant to m x . 

In Fig. [5j we also map direct detection limits into the mz'-m x plane, assuming g q = 
g x = 0.5. For the FV and SV models, direct detection is much powerful than collider 
searches except for m x < 10 GeV. Hence the direct detection and collider detection can 
be complementary to each other in the parameter space. For the FA model, the collider 
detection will significantly improve the limits from the direct detection. 

Since the FA model involves a massive vector boson coupling to non-conserved axial vector 
currents, it has a dangerous UV behavior. When rriz' tends to zero, the DM production cross 
section is essentially proportional to g q g^m q m^/(m A Z iQ 2 ) and would blow up. Consequently, 
as shown in Fig. [5^b), the sensitivity for g q = g x > 0.3 is unusually improved for rriz' < 
50 GeV. In fact, the simplified model falls down in this case. The problem could be solved 
in a UV complete model, where more degrees of freedom are introduced and mz> is no 
longer an arbitrary parameter but given by some physical mechanisms, such as spontaneous 
symmetry breaking. As mz> tends to zero, the DM production process in the FA model may 
violate the unitarity condition. In order to estimate the invalid region of this model, we 
derive a unitarity bound from the parton-levcl process qq —> XX with an appropriate center- 
of-mass energy of the quark pair. The detailed derivation is described in Appendix O For 
g q — g x — 1, the unitary bound for bb —> XX are shown in Fig. |5|T) by a region with light red 
color. In this region, the FA model cannot be a correct description and the collider limits 
derived in this approach are meaningless. 

Finally, let us carefully inspect the connection between the collider searches and the DM 
relic density. The observed DM relic density recently given by the Planck collaboration 
is Xl x k 2 = 0.1173 ± 0.0031 59]. Assuming DM is thermally produced, the regions with 
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(a) FV model 


(b) FA model 



(c) SV model 

FIG. 6: Estimated 90% C.L. limits (dashed lines) on the FV (a), FA (b), and SV (c) models in the 
mz'-m x plane for the mono jet + channel at a pp collider with y/s = 50 TeV, in comparison with 
the regions allowed by the observed DM thermal relic density. Red/green/blue lines correspond to 
g q = g x = 0.1/0.3/0.5. Regions enclosed by solid curves predict a relic density smaller than the 


Planck+WMAP observation Q Y h 2 = 0.1173 


59]. 


small annihilation cross section may overproduce DM in the early Universe and thus be 
excluded. On the other hand, the regions with DM underproduction may still survive in 
some situations, for instance, there exist other DM components, or DM particles are non- 
thermally produced from decays of new heavier particles. 

We show in Fig. [6] the regions where the estimated thermal DM relic density is smaller 
than the observed value. In the FV model, the monojet search at a pp collider with yfs = 
50 TeV is expected to explore most of the allowed regions for g q = g x > 0.3. In the FA 
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model, some allowed regions with m x > m t cannot be tested by the collider search. When 
estimating the relic density in the SV model, we also include the annihilation contribution 
from the interaction term g x X*xZ l ll Z' ll . I which is naturally induced in the U (1) gauge invariant 
extension of the SV model. This contribution can be significant for m x > mz>■ For g q = g x = 
0.5, the monojet search can cover the entire region allowed by the relic density observation 
in the SV model. 

If there are some other interactions that are not included in our simplified models, the 
thermal DM annihilation can be enhanced and reduce the overproduction of DM particles. 
In this case, the monojet search may not be affected by the additional interactions, but 
the regions allowed by the observed relic density would be enlarged. For instance, there 
may exist some Z'-lepton interactions which enhance DM annihilations without affecting 
DM productions at pp colliders. Other possibilities may arise from the UV extension of 
the simplified models. If the mass of Z' and x are obtained via a spontaneously symmetry 
breaking mechanism, an additional Higgs boson h' may exist and mediate a new annihilation 
process xx!XX* ~ ► h' —> Z'Z'. In this case, the correct relic density may be more easily 
achieved. 

V. CONCLUSIONS AND DISCUSSIONS 

In this work, we consider a class of DM simplified models and investigate the prospect 
of future pp colliders. These models contain a new vector boson connecting DM particles to 
quarks via vector or axial vector interactions. In particular, we study the monojet channel 
at future colliders with collision energies of 33, 50, and 100 TeV. 

In the simplified models, full kinematics and topologies can be well investigated. The 
DM pair-production cross section depends on the DM particle mass, the mediator mass 
and the mediator couplings. Future colliders have capabilities to search for the mediators 
with masses of d(TeV). If m x < mz'/Z, the mediator is on-shell produced and then decays 
into DM particles. In this case, the DM production rate is resonantly enhanced. On the 
other hand, if m x > mz'/2 the sensitivity would drop quickly since the mediator is off-shell 
produced. As a result, future colliders may be more sensitive to a heavier on-shell mediator 
than a lighter off-shell mediator. 

We also compare our expectation of the collider sensitivity with constraints from the 
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direct and indirect DM detection. Although the collider phenomenologies of vector and axial 
vector interactions are similar, signal features in other DM detection experiments can be 
quite different. For the vector interaction, collider searches have stronger capability to detect 
light DM with masses smaller than ~ 10 GeV. On the other hand, the collider sensitivity to 
the axial vector interaction can be better than direct searches by several orders of magnitude. 
Moreover, collider searches would be much more sensitive than indirect searches for DM 
particles with masses smaller than the mediator mass. Furthermore, future collider searches 
could explore the bulk of the parameter space that is allowed by the observed DM relic 
density. 

Future pp colliders with very high collision energies will be sensitive to very heavy particles 
in new physics models. In general, other energetic objects, such as charged leptons, photons, 
and 6-jets, can be used to trigger DM signatures. With a very high luminosity, one can 
even utilize reconstructed gauge bosons, Higgs bosons, and top quarks to search for DM 
particles. These kinds of searches may benefit from small SM backgrounds at a high energy 
scale. The challenge is how to accurately measure objects with large px of ~ C?(TeV) 
and how to efficiently reconstruct particles that are highly boosted. The DM study on these 
collider signatures will be a great complement to the monojet search and other DM detection 
experiments. 
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Appendix A: DM-nucleus scattering cross sections 

Direct detection experiments measure the nuclear recoil energy induced by the DM- 
nucleus scattering, where the momentum transfer is far smaller than the reduced mass of 
the DM-nucleus system. In the limit of zero momentum transfer, the mediator in the DM- 
quark scattering can be safely integrated out. We can use effective Lagrangians to describe 
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DM-quark interactions as 

£ fv = -^2G xq X'fx<hnQ, ( Al ) 

g 

= ~^G xq x ( A2 ) 

g 

Cfy - - V G x ,{x'd“ x - (d"x’)x}<h„q, (A3) 

q 

where G xq = g q g x /m 2 z , is the effective DM-quark coupling. The vector interactions (lAlj) and 
(jMD lead to SI scatterings, while the axial vector interaction (lA 2 j) leads to SD scatterings. 

DM-quark interactions induce DM-nucleus interactions. The DM-nucleus scattering cross 
sections a xA are given by 
2 2 

<tsi. x a = —— x — —— [ZGy, xp + (A — Z)G \ Xxn ] 2 for the FV and SV models, (A4) 
7 T v n \ A 171 a ) 2 

ITT/ 2 o 

^"sd ,\A = —7 - 7 - - (SpG A , XP + S„G A>xn ) for the FA model. (A5) 

7r \ m x ± m A) 

Here m A is the nucleus mass. Z and A are the charge number and mass number of the 
nucleus, respectively. S£ ( S£) is the expectation value of the spin contributed by protons 
(neutrons) in the nucleus. The effective DM-proton and DM-neutron couplings G xp and G xn 
are given by 


Gy, xv ~ 2G V „ + G 

G A, xv 


G\ 


— G yu + 2 Gh 


q=u,d,s 


G x , AJ. Gx, X n — E G xg AJ, 

q=u,d,s 


(A6) 

(A7) 


where the nucleon form factors for the axial vector interaction are A d = A£ = 0.842 ±0.012, 
A™ = A p d = -0.427 ± 0.013, and A” = A= -0.085 ± 0.018 60], 


Appendix B: DM annihilation cross sections and relic density 


DM annihilations could induce high energy cosmic rays, gamma rays, and neutrinos, 
which may be observed by indirect detection experiments. The DM annihilation cross sec¬ 
tions into quarks in the simplified models are given by 


^FV.ann (XX q<l) = 


Pq 


l(9q9xY 


127r/3 x (s — m|,) 2 ± 


2 m 2 

1 +- 1 


2 m 2 

1 +-* 


(Bl) 
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0FA,i 


i(xx ^qq) = 
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Pxx* ^qq) = 


P q P x 
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2m 2 
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m 
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(B2) 

(B3) 


where s is the squared center-of-mass energy of a DM particle pair, /3 X = 1 — 4m 2 /s, 

and P q = 1 — 4m 2 /s. DM particles can also annihilate into Z' pairs via f-channel and 

?j-channel if m x > For the SV model, we also consider the contribution from the 

interaction term 5 , x X*xZ^ t Z /At introduced in a U(l) gauge extension of the simplified model. 

Taking into account the velocity distribution of DM particles, the quantity directly con¬ 
necting to indirect searches is the thermally averaged annihilation cross section ( cr ann VM), 
where vm = a/ {pi ■ P 2) 2 — m 2 m|/ (EiE 2 ) is the Mpllcr velocity. As pointed in Ref. 6 ll |. 
instead of directly calculating {cr a . nT , Vm), one can conveniently calculate (cr ann t’ re i) in the lab¬ 
oratory frame and give the same result. Here the laboratory frame means one of the two 
initial particles is at rest, and ry e i is the relative velocity between them. 

In the laboratory frame with the low DM velocity, s can be expanded as s = 4m 2 + 
m 2 v 2 + |m 2 u 4 + 0{v 6 ), where v = u re i = P x { 1 — 2m 2 /s) _1 . Substituting this expression into 
(T ann , we can expand cr ann w as a + 6u 2 + (2(u 4 ). For a Maxwell-Boltzmann velocity distribution 
parametrized by temperature T, we have (a ann 7>) = a + 6 bx~ l + 0{x~ 2 ) 1 where x = m x /T. 
We list below the coefficients a and b in the simplified models. In the FV model, 
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where the step function means the annihilation channel 
In the FA model, 

g x (l-m 2 z ,/m 2 ) 3 / 2 

“ fa = - mz,) 


XX Z'Z' opens only if m x > my. 
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a fa for DM annihilations into quarks are proportional to m 2 , implying that the annihilations 
are hclicity suppressed in s-wave. In the SV model, 


g^l-m 2 z ,/m 
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Z'J 


DM annihilations into quarks are of p -wave in the leading order, because a pair of scalar 
particles cannot form a vector state without orbital angular momentum. 

The evolution of thermal DM density can be determined by the Boltzmann equation 

dn v 


dt 


+ 3 Hn x = - (a ann u) (n x n x - n® q n| q ) = - (a ann v) [(n x ) 2 - (n® q ) 2 ], 


(BIO) 


where H = ^87t p/ (3M 2 j) is the Hubble rate, p is the energy density in the Universe, and 
M p i is the Planck mass. n x (n x ) is the number density of the DM particle (antiparticle). 
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The superscript “eq” represents the thermal equilibrium. Here we do not consider DM 
particle-antiparticle asymmetry and assume n x = n x . 

In principle, one can numerically solve the Boltzmann equation and get the DM relic 


density. Here we use an approximate method to obtain 62, 63] 

3 


nji 2 = 2 x 1.04 x 10 9 GeV 


-i 


T n 


x f 


2.725 Kj M pl ^gjxf)(a + 3 b/xf )’ 


(Bll) 


where Xf = m x /Tf ~ (9(10), Tf is the DM freeze-out temperature, T 0 = 2.725±0.002 K 64] 
is the present CMB temperature, and g*(xf) is the effectively relativistic degrees of freedom 
at the freeze-out epoch [65]. 


Appendix C: Unitarity bound 


The unitarity condition for a 2-body inelastic process can be expressed as 

1 


\af e \s) 


< 


(Cl) 


2 VKs 

where /3(s, in) = y/l — 4m 2 /s, m in (m out ) is the mass of either of the two incoming (outgoing) 
particles, and s is the center-of-mass energy of the system. a) nel (s) is the j-th partial wave 
coefficient of the invariant amplitude A4 for the process: 
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a? el (s) = 


dd sin 6Pj (cos 9)M mel (s, cos 9 ), 


J 32t r 
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where Pj(x) is the j'-th Legendre polynomial. 

In the FA model, the 0-th partial wave coefficient for the process q~q~ —> X-X- (th e 
minus sign means a hclicity of —1) is given by 


.inel / 
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9q9x 


m q m x 
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I'Z’ Iff ouuz'J- Z' \ 

Apparently, this process violates the unitarity condition when rriz' \Z§. 

In order to obtain the unitarity bound on DM particle pair productions at pp colliers, 
one challenge is how to determine s, which varies event by event and is affected by parton 
distribution functions and associated jets. There are several estimates of .s used in the 
literature, such as the center-of-mass energy of the collider { 17 } ]. yj 4 m x + pj, Q S], and 


(C3) 


an estimate of the invariant mass of the DM particle pair in some fraction of events 


66]. In 
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this work, we adopt the strategy in Ref. 19] and take s as the averaged momentum transfer 
weighting with parton distribution functions: 


„ _ ,Q2\ _ I dx^jfqjx^fgjxz) + fq{X2)fq{Xl)]&(Q ~ 2 TTl^Q 2 
^qj dxidx 2 [f q (x 1 )fq(x 2 ) +f q {x 2 )fq{x 1 )]Q(Q-2m x ) 

where Q 2 = (p q + Pq — Pj) 2 = X\X 2 s — y/sp'r(xie~' n + x 2 e v ) is the squared momentum transfer 
to the DM particle pair in the process qq — > XX3- Pt and r) are the transverse momentum 
and pseudo-rapidity of the jet j, respectively. The unitarity bound demonstrated in Fig. [5] 
is derived from the process —> X-X-- We have set rj = 0 and px as the cut threshold 

in order to give a conservative bound. 
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